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Abstract 

We  present  a  new  method  for  approximating  the  partition  function  of  2D  Ising  model.1'  usin^ 
a  transfer  matrix  of  order  2".  For  //  =  40  our  current  program  tock  about  20  seconds  on 
a  Sparc  station  to  obtain  4  correct  decimals  in  the  lop  two  eigenvalues  and  •">  minutes  for 
0  correct  decimals.  Eigenvectors  were  computed  at  ihe  same  time.  The  lemperature  was 
within  f'/T  of  critical. 

The  main  idea  is  to  force  certain  entries  in  veclors  to  have  the  same  values  and  to  find 
the  crudest  representation  of  this  type  that  delivers  the  required  accuracy.  At  no  time  doe^ 
our  program  work  with  vectors  with  2"  entries. 


1  Introduction 


The  Ising  model  was  proposed  to  explain  properties  of  ferromagnets  hut  since  then  it  has 
found  application  to  topics  in  Chemistry  and  Biology  as  well  as  in  Physics.  For  any  reader 
unfamiliar  with  the  mode!  an  excellent  introduction  targeted  at  a  genera]  audience  is  [C'ip*7]. 
The  remainder  of  this  section  assumes  some  knowledge  of  the  so  called  transfer  matrix.  This 
paper  presents  a  numerical  method  for  computing  properties  of  the  2D  Ising  model  for  given 
parameter  values  such  as  magnetic  field  strength  B.  temperature  T  and  coupling  constants 
.7. 

There  are  two  avenues  leading  to  such  calculations:  combinatorial  and  algebraic.  Our 
method  is  in  the  second  category  which  makes  use  of  a  transfer  matrix  M„  associated  with 
a  semi-infinite  helical  grid  of  "spins"  or  "sites"  with  n  of  them  on  each  circular  band.  One 
form  of  M„  for  u  =  3  and  n  =  4.  with  the  field  strength  B  normalized  with  respect  lo  t he 
coupling  constant  ■)  is  as  follows: 


1 


where  (with  appropriate  normalizations) 


«  =  b  =  <~Bir  and  r  = 

Tlie  attractive  property  of  .17,,  is  that  it  is  a  nonnegative  irreducible  matrix  whose 
dominant  eigenvalue  (called  the  Perron  root)  is  the  wanted  partition  function  per  spin. 
Thus  it  is  only  necessary  to  approximate  this  eigenvalue  to  the  desired  accuracv  although 
the  associated  eigenvectors  are  also  useful  in  approximating  quantities  of  physical  interest. 
Moreover  .!/„  is  exceedingly  sparse:  it  has  exactly  2  non-zero  entries  per  row  (and  column) 
arranged  in  a  regular  pattern.  There  is  only  one  difficulty:  .17,,  is  of  order  2”  and  we  are 
interested  in  the  case  n  —  We  know  of  no  calculations  with  u  >  20  up  till  now. 

Our  approach  uses  a  finite  family  of  orthogonal  md/rinl  vectors,  and  approx¬ 

imates  the  lop  two  column  and  row  eigenvectors  of  .17,,  from  the  subspace  spanned  hv 

S„.i. 

Step  0.  Init ialize  /  to  1 . 

Step  1.  Represent  in  compact  iorm.  the  orthogonal  projection  I'  of  the  transfer  matrix 
.17,  onto  the  subspace  span( 

Step  2.  (  ompute  the  two  largest  eigenvalues  and  the  associated  column  eigenvectors  of 
I’.  1  hese  are.  in  a  sense,  the  best  approximations  from  the  given  judicial  subspace 
spain«S„./).  However  they  may  noi  be  good  enough. 

Step  3.  Lvaluate  residual  norms,  condition  numbers  and  associated  error  bounds  and  es¬ 
timates.  If  the  estimates  are  satisfactory  then  compute  the  required  properties  of  the 
model  and  stop.  Otherwise  return  to  Step  ]  with  the  ne.xi  member  of  each  family,  i.e. 
increase  /  by  1. 

Our  goal  is  to  creep  up  to  the  coarsest  of  our  vector  representations  that  permits  ap¬ 
proximations  of  the  desired  accuracy.  This  minimal  representation,  which  is  noi  known  in 
advance,  gave  in  the  name  for  our  approach. 

Note  that  the  difficulty  lies  not  in  .17,,  itself  but  in  the  the  representation  of  vectors 
in  R-  .  Indeed  the  special  structure  of  ,17r,  would  permit  evaluation  of  .!/„?•  for  any  2r- 
dimensiona!  veclor  r  with  great  efficiency.  However  a  procedure  that  cosis  0(2")  may  be 
too  much  when  »  is  large  and  our  central  problem  is  the  representation  of  vectors  in  R-”  . 

Sparse  vectors  occur  in  sparse  matrix  work  and  N.  Fuchs  [FucXOj.  when  applying  the 
Power  Method  to  .!7„.  keeps  only  the  largest  1000  eniries  of  each  vector.  This  device  is 
satisfactory  deep  within  ihe  ferromagnetic  region  of  the  model.  However  after  studying 
the  Perron  vector  in  cases  near  the  critical  temperature  we  found  that  it  contained  almost 
no  small  entries.  In  different  language,  every  configuration  in  the  "spin”  array  contributes 
significantly  to  the  partition  function. 


A*  a  substitute  for  sparsity  we  propose  lo  limit  the  number  of  distinct  values  that  ran 
occur  among,  a  vector's  components.  YVe  do  this  by  means  of  a  family  of  "indicial  veclors". 
Here  is  a  sketch  of  the  idea.  More  details  are  given  in  Section  2. 

A  vector  in  RJ  may  be  thought  of  as  a  function  on  {1.2 . 2"}.  What  we  call  an 

indicial  function  is  really  a  partition  of  this  index  set  into  disjoint  subsets  on  each  of  which 
the  vector  is  constant.  Thus  the  vector  takes  on  fewer  than  2’’  distinct  values,  perhaps  only 
a  few  million  of  them.  This  sort  of  vector  recalls  H.  Lebesgue  s  approach  to  integration 
via  step  functions.  For  a  given  partition  /  the  set  of  all  representable  veclors  forms  a 
subspace  Sj  of  RJ  .  We  bow  to  the  influence  of  computer  science  and  start  counting  at 

0.  If  {c0 . denotes  the  standard  basis  and  if  { 1  j.  93. 214.  NGb}  is  one  subset  in  the 

partition  /  then  *  j=,  ■+-  c93  +  c_.]4  +  is  one  member  of  a  natural  orthogonal  basis  for  Sr. 
In  other  words,  the  natural  basis  vectors  of  R-’  are  aggregated  according  to  /  to  produce 
an  orthogonal  basis  of  S/.  An  important  feature  of  our  approach  is  that  these  basis  vectors 
ate  never  represented  explicitly  in  the  computer.  Careful  index  manipulation  takes  their 
place.  Moreover  our  choice  of  /  yields  a  manageable  representation  of  the  projection  /’<•  of 
M„  onto  Sr.  Pj  is  nonnegative  and  irreducible.  Pi  is  not  as  sparse  as  M„  but  we  hold  it 
in  a  compact  form  that  permits  the  efficient  formation  of  P<>r  for  appropriate  ir. 

There  is  some  freedom  in  the  choice  of  the  family  of  /' s.  Our  /' s  are  a  compromise 
between  physics  and  the  very  special  structure  of  M„.  Details  are  given  in  Seciion  2. 

The  next  task  is  to  find  the  Perron  vectors  of  Pr.  Recall  that  the  lop  two  eigenvalues 
of  .1/,,  coalesce  as  ihe  temperature  becomes  critical.  We  have  used  two  approaches; 

(ala  block  power  method  with  a  block  size  of  2. 

lb)  a  nonsvmmetric  Lanczo-  code. 

The  details  are  given  in  { PH f H ] .  It  turns  out  that  it  pays  to  compute  the  two  largest 
eigenvalues  together  with  their  column  and  row  eigenvectors.  The  reason  that  conven¬ 
tional  techniques  such  as  these  are  appropriate  is  that  with  our  current  indicial  functions 
f.  dim  Sj  =  ()[ii~2~])  and  so  P>  is  of  modest  order.  In  addition  we  form  and  compuie 
similar  quantities  for  ();•■  the  (orthogonal  l  projection  of  onto  an  associated  subspace 
Sj-.  The  extra  information  from  Q j<  allows  us  to  compute  an  approximate  Perron  row 
vector  y"  to  match  the  Perron  column  vector  .r  for  Pj.  Pj  and  Qp  share  the  same  Perron 
root.  Fortunately  Qp  is  diagonally  similar  to  Pj  and  need  not  be  represented  explicitly. 

We  would  prefer  to  use  the  oblique  projection  of  .\/n  onto  the  pair  of  subspaces  {Sj.Sp\ 
but  we  have  not  yet  found  a  convenient  (sparse)  representation  because  some  of  the  canonical 
angles  between  Sj  and  Sp  equal  tr/2  and  this  fact  complicates  the  representation. 

Associated  with  the  vectors  x  (Pjx  -  x~)  and  ym  ( y'Qp  =  ~y‘ )  are  vectors  zj  £  Sj 
and  w’j  €  Sp  that  approximate  the  eigenvectors  we  seek.  It  is  essential  to  be  able  to  bound 
or  estimate  the  accuracy  of  our  approximate  eigentriple  »•}). 

Fortunately  by  using  our  special  bases  in  Sj  and  Sj'  appropriately  we  can  compute 


(exactly  in  exact  arithmetic)  the  associated  residual  vectors 

rj  :=  M„zj  --/'/•  ■'•}'■=  trf  ~  u-j'-j' 

and 

-7  :=  «■/-//( ll,r/ll-l!~/ll- )■ 

Although  ij  t  R-”  .  .«/.  €  R*  we  can  accumulate  |jry||*  and  II--/* ||“  and  irj  during  the 
computation  of  zj  and  and  thus  avoid  ever  having  to  store  them.  This  is  a  key  feature 
of  the  efficiency  of  our  method.  From  ||r||.  |!-s||.  and  —7  we  can  compute  error  bounds  and 
error  estimates.  This  is  discussed  in  [PH9lj. 

It  is  like! v  that  our  error  estimates  indicate  that  zr.  irj>  and  ~j  are  not  sufficiently  accu¬ 
rate.  In  that  case  we  pick  the  next  indicial  function  f  in  our  family  so  that  /  is  a  refinement 
of  f  and  St  C  Si.  dim  Si  -  2  dim  Sj.  Then  we  repeat  the  cycle  of  approximations  until 
the  accuracv  requirement  is  met  or  our  resources  are  exhausted.  This  is  noi  an  iterative 
method  because,  in  a  finite  number  of  steps,  ilie  indicial  function  becomes  the  identity. 

I5y  creeping  up  to  adequate  approximations  from  below  we  ensure  that  we  end  up  wiih 
the  coarsesi  indicial  f  net  ion  that  meets  the  given  tolerance.  In  this  way  do  we  achieve  the 
minimal  representation,  from  our  family,  that  gives  our  method  it'  name. 

It  is  worth  repeating  t ha  1  at  no  time  in  the  cycle  do  we  need  to  store  a  vector  with  2 
components. 

Quantities  of  interest  are  usually  partial  derivatives  of  the  partition  function.  If  we  used 
differences  to  estimate  derivative'  that  would  sharply  increase  the  required  accuracy  oi  our 
approximations.  Fortunately  S.  Gartenhau'  [Garvy  and  X.  Fuchs  fFucSbj  have  shown  that 
some  of  the  quantities  of  interest  may  be  expressed  in  terms  of  r  and  it ■*  and  so  there  is  no 
need  to  use  differences.  This  increase'  the  scope  of  our  approach  significant ly. 


A 


2  Projections  of  duodiagonal  matrices 


Figures  4  and  -1  show  that  the  Perron  eigenvector  of  may  he  approximated  quite 
well  by  vectors  in  which  certain  positions  are  constrained  to  carry  the  same  value.  The 
challenge  is.  of  course,  to  specify  in  general  ihe  right  positions.  Our  approach  makes  heavy 

use  of  the  binary  representation  of  the  numbers  0.1 . 2"  -  1.  In  particular  we  always 

use  //-hits  in  a  representation.  Thus  the  positions  in  a  2'*- vector  are  indexed  by  the  bit 

strings  {00  •  •  -0.  00  •  •  -01 . 1 1  •  -  •  1 1 }.  If  we  think  of  these  strings  as  configurations  of  a 

ID  Ising  model  then  we  obtain  equal  contributions  to  the  (total)  partition  function  from 
all  configurations  (strings)  which  have  the  same  /,■  (the  number  of  l‘s)  and  1  (the  number 
of  transitions).  Thus  we  can  group  the  bit  strings  by  the  values  of  (k.1).  For  //  =  -1  this 
partition  is  shown  in  Figure  1. 
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Figure  1:  Indicia!  sets  and  basis  vectors  for  bit  strings  of  length  //  =  4. 


The  approximation  is  not  bad  at  the  visual  level  bul  may  wen  not  be  adequate.  We 
need  a  systematic  way  to  refine  the  partition  induced  by  k  and  t.  Our  choice  may  not  be 
optimal  but  it  has  the  practical  virtue  of  exploiting  the  duodiagonal  structure  of  M„.  We 

use  the  last  /  bits  in  the  string,  for  /  =  1.2 . as  a  refining  parameter.  For  large  enough 

/  the  original  index  set  is  recovered.  The  partition  for  n  =  o.  /  =  1  is  shown  in  Figure  2. 
In  each  figure  the  last  column  lists  a  set  of  orthogonal  vectors  in  R-’  which  we  call  S,,j: 
v  =  4.  /  =  0  for  Figure  ] .  //  —  5.  /  =  1  for  Figure  2.  A  careful  analysis  (see  [Hen91])  shows 
the  cardinality  of  this  set: 


kV/i  =  2'(1  + 


•V/i  ,  (»  -  1+  -  1),  ^  „2.V-i 


<  7/ *2 


Our  method's  utility  depends  on  obtaining  adequate  accuracy  for  small  values  of  /.  However 
defining  S„j  is  not  enough.  We  need  to  derive  the  projection  P^ ;  of  Mn  onto  span(«Sti./l 
without  actually  using  the  indicial  vectors  exhibited  in  the  figures  above.  Instead,  by 
analysis,  we  determine  a  priori  the  positions  and  values  of  the  nonzeros  in  P^,. 

The  duodiagonal  form  of  Mn  is  essential  to  the  analysis.  The  key  fact  is  that  the  action 
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Figure  2:  Indicial  sets  and  basis  vectors  for  u  =  5  and  /  =  1. 
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of  M„  on  any  vert  or  in  S„j  ran  be  expressed  as  a  linear  combination  of  either  2  or  1  vector' 

in  $n.i+ 1-  To  illustrate  the  idea  we  show  how  to  obtain  a  column  of  Pft  t  for  »  —  ■’>.  /  =  1. 

Let  w-  denote  the  last  I  bits  of  an  n-string  and  let  denote  the  indicial  vector  specified 

by  A.  and  /.  e.g.  {00010.00100.01000}  is  the  indicial  set  with  A  =  1.  1  =  2.  ~  =  0  (with 
n  =  .)./=  1).  Thus  .tu  i.i  =  f  i  +  +  eg.  The  duodiagonal  structure  requires  that  the 
uonzeros  in  column  j  occur  in  rows  2 j  mod  2'  and  (2 j  +  1 )  mod  2". 

•'W/I't' 0.1.2  =  d/ C  j  Mn  C  _j  -f  Mr.  <  8 

=  (  (I (  4  +  bf  S  )  +  (  etc  S  +  hf  ;)  )  +  (  AlC  1  < .  +  CC  1  7  ) 

=  (l{‘  A  +  f  £  )  +  6(  f  .s  +  f  9  )  +  H-  +  1  7 

=  w-t‘oo.i.2  +  ot.2.3  +  Au‘oo.i.1  +  c.i'01. 2.2- 

See  the  matrices  illustrated  in  the  Introduction  for  the  meaning  of  «.  b.  and  c.  Next  we 
must  determine  the  triples  such  thal  the  inner  ]>iod u ct 

does  not  vanish.  We  return  to  our  example  above  and  read  off  tlie  results  because  all  indicial 
vectors  are  pairwise  orthogonal. 


( -too. i.i-  0.1.2)  =  00.1 . j ! I ”  =  2(i 

{•cm. 1.2-  -W„.r0.i.j)  =  A»||j'oi.i.j||*  =  2b 

b'ou.i.i- -M,. -'0.1.2)  =  ^ll-'oo.i.i  li'  -  b 

(■I'M. 2.2-  -1/,. -'0.1.2')  =  ^il-t'oi .2.2 =  ( 


All  other  entries  in  column  (0.1.2)  vanish.  Since  the  indicial  vectors  are  not  normalized  a 
diagonal  scaling  is  necessary.  The  result  is  that  the  -1  nonzeros  in  column  (0.  1.2)  of  , 
are  just  a.b.b.c  and  they  occur  in  the  indicated  rows.  Thus  P  is  easily  siored.  in  compact 
form,  in  a  rectangular  array  as  discussed  in  the  next  seel  ion.  Some  other  columns  of  P  will 
have  only  2  nonzeros.  In  the  general  case  the  precise  form  of  P ; ;  is  not  easy  to  determine. 
However  the  analysis  has  been  done  and  the  details  may  be  found  in  [Henbl],  The  result 
is  that  P  may  be  computed  with  arithmetic  effort  proportional  to  |<Sn./|  and  with  storage 
proportional  to  i5„./+i|.  By  means  of  this  combinatorial  analysis  the  largest  two  eigenvalues 
of  M„  may  be  approximated  by  the  largest  two  eigenvalues  of  /'{},  and  this  can  be  tackled 
by  more  conventional  techniques. 

Here  is  a  sample  of  the  sizes  of  the  sets  and  the  maximal  number  of  ones  in  an 
indicial  vector  x  = 
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Table  1:  Combinatorial  properties  of  suffix-based  indicia]  sets 
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3  Implementation  issues 


We  discuss  the  data  si  induces  used  for  indicia!  sets  and  indicia!  vectors  ./•_;.  as  well  a- 
efficient  algorithms  for  manipulating  coefficient  vectors. 

3.1  Representation  of  vectors 

In  order  to  turn  t lie  set  S„.i  of  indicia!  vectors  r  into  a  basis  a  total  ordering  on  tin* 

index  triples  (w\  A-.f)  with  |— j  =  /  is  needed.  Recall  that  — •  is  the  /-bit  suffix  of  an  //-string. 
Such  a  triple  is  hyiliniaU  if  the  corresponding  index  set  is  nonempty  (or.  equivalently,  if 
J -.h.t  is  a  vector  in  S„j).  The  ordering  on  legitimate  triples  uses  the  number  ri _•  |  of  which 
~  is  the  binary  representation: 

(--'.a •'./') 

if  and  only  if. 
i  a  )  /•(_'  )<  i.  or 


(b)  /•(-')  =  and  A-  <  A',  or 

(c)  /•(_•)  =  /•(-■').  A-  =  1/  and  1  <  1' . 

We  use  for  the  ranking  of  i  under  ibis  ordering,  i.e.  <J>  maps  („•./,•./)  into  a 

nonnegative  integer. 

Let  .V  denote  the  matrix  whose  column-  are  the  vector-  in  S,.j  in  the  prescribed  order. 
Then  any  vector  g  in  span(5,,.;)  may  be  written  as 

<1  =  -Xg 

where  g  is  the  suffix-bast  (I  cot  ffirif  ul  rtclot  of  g.  Note  that  g  ~  R*  while  g  ■£ 

It  is  tempt  in  g  to  use  1  as  an  index  so  that  any  g  could  be  held  in  a  31)  array  of  real 

numbers.  The  trouble  with  this  scheme  is  that  there  are  holes  i illegitimate  triples)  and  so  it 
is  a  bit  wasteful  of  storage  and.  even  worse,  every  access  to  the  array  must  check  whether  the 

index  is  legitimate.  Instead  we  simply  map  the  legitimate  triples  onto  {0.1 . -  1} 

using  $. 

The  first  task  is  to  obtain  4>  given  (-.A './).  The  definition  of  -</  shows  that  1  plays  a 
minor  role  and  so  we  only  need  a  2D  array,  based  on  r(-)  and  k.  which  we  call  nuhx.  The 
needed  entries  in  the  21  >.  I »  +  1 )  array  index  may  be  computed  initially.  Then 

4» ( _' .  A- .  / )  =  /nc/c.r[r(w')][A']  -f  1. 

The  next  task  is  to  represent  the  inverse  mapping:  given  a  value  for  $  find  k.  and  i. 
This  is  easily  accomplished  with  three  ID  arrays  of  length  |«S„,i|.  one  for  each  of-'.  A\  and 
t.  These  are  initialized  before  the  calculations  begin. 
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3.2  Advantage  of  the  indexing  scheme 


VYlien  computing  t lie  dominant  eigenvectors  of  a  projection  matrix  t lie  array  <j  i>  treaied  a~ 
a  conventional  vector  in  hut  when  we  need  to  know  an  entry  of  g  —  Xg  in  RJ  we 

proceed  as  follows.  To  find  <i(i)  first  obtain  ft.  t lie  n-bii  binary  representation  of  i.  together 
with  it'  /- l)it  suffix  _•  (the  last  /  bits  of  //).  Next  compute  k  and  1  from  //.  Finally  look  up 
4>(w .  k.  1 )  as  described  above.  Then 


</(i)  =  </(4>U-.A.  /)). 

It  is  these  simple  bit  manipulations  and  table  look-ups  that  enable  us  to  avoid  the  use  of 
an  array  with  2"  entries. 

3.3  The  Projection  Matrix 

As  mentioned  earlier  P^  •  has  at  most  -1  nonzeros  per  row  and  so  coniains  somewhat  less 
than  -1; $„/  nonzeros.  We  represent  P  a-  a  sequence  of  packed  columns  in  a  II)  array 
loLproj.  In  addition  two  inn  ex  arrays  ml  and  row  are  needed:  the  /th  nonzero  is  in 
position  i /•»«•(/).  ro/( /))  and  ha*  the  value  ml./noji  /’).  The  vector  P\  r  is  formed  by  taking 
the  linear  combination  of  P  -  columns  with  coefficient  -  given  by  r. 

3.4  Application  of  the  transfer  matrix  to  approximate  eigenvectors 

Although  the  use  of  suffixes  _■  to  refine  our  model  has  no  justification  from  physic*,  ii  doe- 
have  the  virtue  that  the  action  of  M„  on  vectors  in  span! S,,.: )  can  lie  computed  exactly 
(modulo  round  off  error).  Consequently  the  norms  of  residual  vectors  may  be  computed 
without  storing  arrays  of  length  2".  The  reason,  briefly,  is  that  ,\/r,  maps  span(5M.;l  into 
spanu<\,./+i  ••  Details  are  given  in  ( P H 0 J ] . 

3.5  Extracting  Information  from  the  Projections 

Even  though  \$n.l\  <  2”  it  is  essential  to  use  a  fairly  efficient  method  to  compule  the 
two  dominant  eigenvectors  of  P(u  !  and  the  two  dominant  row  eigenvectors  of  (which  i> 
diagonally  similar  to  Fy  ,  and  so  does  not  need  separate  treatment ).  Although  we  only  want 
the  dominant  eigenvector  we  consider  the  calculation  of  two  eigenvectors  to  be  essential  for 
efficiency  when  the  temperature  in  the  Ising  model  is  near  critical  and  there  is  less  and  less 
separation  between  the  two  largest  eigenvalues. 

We  have  tried  two  methods:  (a)  the  block  power  method  with  block  size  2  (called 
subspace  iteration  by  numerical  analysts  and  structural  engineers),  (b)  the  unsymmetric 
Lanczos  algorithm. 

A  block  Lanczos  algorithm  with  block  size  2  would  probably  be  more  efficient  but  we 
have  not  developed  a  code  for  that  yet.  Indeed  the  unsymmetric  Lanczos  code  is  not  yet 
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a  standard  method  IhH  in  our  experiments  it  heroines  increasingly  better  than  the  h]ork 
power  method  as  n  increases  beyond  15. 

More  details  about  the  implementation  and  our  error  estimates  are  given  in  (PHOij. 
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4  Numerical  results 


Here  are  the  results  from  a  preliminary  code  using  the  nonsy linnet ric  Lanczos  algorithm. 
For  the  hardest  case,  v  -  .'10  and  temperature  within  3*7  of  critical,  it  took  about  20  seconds 
on  a  Sparc  station  to  obtain  the  partition  function  to  3  decimal  digits,  and  about  5  minutes 
to  obtain  5  decimal  digits.  In  the  tables  below.  GRQ  is  the  generalized  Rayleigh  quotient 
y~  M „x / y~ x.  The  temperature  T  —  1.0  is  deep  within  the  ferromagnetic  region.  T  =  2.2  is 
within  3*7  of  the  critical  temperature. 


/ 

approximat  ion 

CRQ 

dim 

time  ( s 

2 

3. 01X9X07014 

(2.7  >  Hr") 

3.51X9x22207 

( —  1  .X  x  10~‘:;) 

14x 

0.x 

3 

3.51X9X42995 

(2.9  >  10"7) 

3.51x9x377x2 

(-2.4  >  10"7) 

232 

1.3 

4 

3.51 X9X39750 

\ 

V 

y 

o 

1 

<■/ 

3.51X9X39519 

(-0.2  >  10"s) 

352 

2.0 

Table  2:  Results  for  »  =  10.  B  =  0.0001.  7  =  1.0  (true  eigenvalue  -  3.51x9x40135) 


/ 

approximat  ion 

(IRQ 

dim 

1  ime  ( s 

•> 

2.5925207940  (2.3  •  K)-'*) 

2.5922407533  (-5.1  > 

10“' 1 

14x 

0.x 

3 

2.5923300340  (4.4  >  10"'] 

2.5921X03040  (-1.1  > 

10”^  ! 

232 

1.5 

4 

2.5922000120  (-2.0  *  10"'') 

2.592220()044  (  -(,.0  . 

10“'  l 

352 

2.4 

Table  3:  Results  for  u  —  10.  B  =  0.0001.  T  -  2.2  (true  eigenvalue  =  2.5922922453 ) 


/ 

approximation 

GRQ 

approximation  -  GRQ 

dim 

time  1  s 

2 

3.51X9X02741 

3.51x9759x7.x 

4.3  >  10"*' 

(jxx 

5.2 

3 

3.51x97x0552 

3.51X9731775 

4.9  y  10-'' 

1232 

x.3 

4 

3.51X9775223 

3.51X9777525 

-2.3  >.  nr7 

2192 

17.0 

.) 

3.51x9770100 

3.51X9775001 

5.0  x  10"^ 

3X72 

35.7 

(1 

3.51X9770241 

3.51X9770145 

9.0  >  10";* 

07x4 

71.7 

l 

3.51X977040X 

.1X97771X4 

-7.x  X  10"s 

11770 

132.0 

Table  4:  Results  for  v  =  20.  B  ~  0.0001.  T  -  1.0 


/ 

approximation 

(IRQ 

approximation  -  (IRQ 

dim 

t  ime  ( ' 

2 

2.5*75104697 

2.5*7301 1057 

2.2  >  10"  4 

0** 

0.* 

3. 

2.5X73559732 

2.5*71*52423 

1.7  >■  10-4 

1232 

5.2 

4 

2.5*72924943 

2.5*72247*** 

0.*  v  10“S 

2192 

11.2 

5 

2.5*6**5053* 

2.5*09010709 

-1.7  >  10"' 

3*72 

17.3 

6 

2.5*7257601* 

2.5*72*09*94 

-2.3  >  10~' 

07*4 

*0.9 

1 

2.5X72473229 

2.3*729*1335 

-5.1  x  10-' 

11770 

70.7 

Table 

•5:  Results  for 

n  =  20.  B  =  0.0001.  T  = 

2.2 

/ 

approximat  ion 

(IRQ 

approximation  -  (IRQ 

dim 

lime  | - 

2 

3.51*979*030 

3.51*9277*29 

5.2  >  10-' 

1 02* 

11.* 

3 

3.51*7421095 

3.51*72710*5 

1.5  x  10"' 

3032 

10.0 

4 

3.51*9705902 

3.51*9734 194 

3.2  >  10-'' 

5032 

50.9 

5 

3.51*9754*09 

3.51*9030*14 

1.2  •  Hr' 

10432 

101.5 

0 

3.51*9707320 

3.51  *9705430 

1.9  -  J  or7 

19204 

21 3.3 

1 

3.31*9774542 

3.51*9775232 

-0.9  >  10"^ 

354  50 

4  72.3 

Table  0:  Results  for 

n  =  30.  B  =  0.0001.  T  = 

1.0 

/ 

approximation 

(IRQ 

approximation  -  GRQ 

dim 

time  (s 

2 

2.5*05*77390 

2.5*04247904 

1.6  >  10"4 

102* 

21.1 

3 

2.5*04495900 

2.5*03307035 

1.1  >  10"4 

3032 

17.5 

4 

2.5*039*93*9 

2.5*03409514 

5.*  x  10~5 

5032 

3*.  3 

5 

2.5*0373*510 

2.5*0342905* 

3.1  x  10"' 

10432 

04.3 

6 

2.5*03033205 

2.5*63020747 

1.2  x  10-6 

19204 

139.1 

i 

2.5*03035130 

2.5*03*31549 

-2.0  x  10-' 

35456 

310.5 

Table 

7:  Results  for 

n  =  30.  B  =  0.0001.  T  = 

2.2 
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Figure  4:  Dominant  column  eigenvector  of  \In.  n  =  8.  B  =  0.0001.  7 
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